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The gravitational aether theory is a modification of general relativity that decouples vacuum 
energy from gravity, and thus can potentially address the cosmological constant problem. The 
classical theory is distinguishable from general relativity only in the presence of relativistic pressure 
(or vorticity) . Since the interior of neutron stars has high pressure and as their mass and radius can 
be measured observationally, they are the perfect laboratory for testing the validity of the aether 
theory. In this paper, we solve the equations of stellar structure for the gravitational aether theory 
and find the predicted mass-radius relation of non-rotating neutron stars using two different realistic 
proposals for the equation of state of nuclear matter. We find that the maximum neutron star mass 
predicted by the aether theory is 12% - 16% less than the maximum mass predicted by general 
relativity assuming these two equations of state. We also show that the effect of aether is similar to 
modifying the equation of state in general relativity. The effective pressure of the neutron star given 
by the aether theory at a fiducial density differs from the values given by the two nuclear equations 
of state to an extent that can be constrained using future gravitational wave observations of neutron 
stars in compact systems. This is a promising way to test the aether theory if further progress is 
made in constraining the equation of state of nuclear matter in densities above the nuclear saturation 
density. 



I. INTRODUCTION 

One of the fundamental challenges of modern physics 
is to solve the cosmological constant problem. This prob- 
lem, in its various forms, has been with us since the 
beginning of the 20th century. At the end of the cen- 
tury, data from Type la supernovae pointed to the fact 
that the expansion of the universe is accelerating [THH]- 
A constant is needed in the general relativistic descrip- 
tion of gravity in order to explain the acceleration of the 
cosmic expansion. One can interpret this "cosmological 
constant" as the energy of the vacuum (popularized as 
dark energy). Quantum field theory predicts a value for 
this energy. However the predicted value is much larger 
than the value found from observational data. This is the 
so-called "old cosmological constant problem" [4 . The 
new problem is why this constant is very small but not 
zero. Another problem is that the values of dark energy 
density and matter density are found to be of the same 
order at the present epoch oddly suggesting that we are 
living in a special time in the history of the universe. 
This is the so-called "coincidence problem" . 

One possible solution to the cosmological constant 
problems was suggested by one of us in [5] , which mod- 
ifies the Einstein equation in the following way so that 
the vacuum does not gravitate: 

(S^r^M = - l -Tg^ + ... , (1) 
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where Q — 4Gjv/3 and Gat is the usual Newton's gravi- 
tational constant. By subtracting the trace of the energy 
momentum tensor on the right-hand side, the Einstein 
equation becomes insensitive to the vacuum energy den- 
sity, Pvac, where 7) t „ = Pvacfl 1 /^ + excitations. As energy 
and momentum are conserved, the divergence of T^„ van- 
ishes. By definition the divergence of G^ v also vanishes 
through the Bianchi identities. Therefore, if we want to 
subtract the trace of the energy momentum tensor, we 
need to add a suitable term to it so that the divergence 
of the right-hand side of Eq. [I] vanishes. It was suggested 
in [5] that this term can be the energy momentum ten- 
sor of a perfect fluid 7^„, which is dubbed "gravitational 
aether". With this term, Eq. [I] takes the form: 

(8irgy 1 G^[g^} = T ^ - \ T 9^ + V> ( 2 ) 

% v = V{U^U v +g liV ), (3) 

where V and are the pressure and four- velocity of the 
gravitational aether. For the right hand side of Eq. [2] to 
be divergenceless, we require: 

It is argued in [5] that the pressure and four-velocity 
of the gravitational aether are dynamically fixed in terms 
of T^ v via Eq. [3] Furthermore, it is consistent with all 
the current precision tests of gravity and cosmological 
observations [5J. 

With the metric being blind to vacuum energy, the 
gravitational aether theory solves the old cosmological 
constant problem. But how does it explain the coinci- 
dence problem? In [7] the authors study static black 
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hole solutions in the gravitational aether and argue that 
aether couples the spacetime metric close to the black 
hole horizon, to the metric at infinity. They then show 
that this can lead to an accelerating cosmological solu- 
tion, far from the horizon. This connection between the 
formation of stellar black holes and the acceleration of 
the expansion of the universe can potentially solve the 
coincidence problem, through Planck-suppressed correc- 
tions in black hole physics. 

All this means that the gravitational aether theory is 
an attractive alternative to general relativity. However, 
the theory is good only if it can make definite predictions 
that are consistent with current and future observational 
results. In astrophysical situations where gravity due to 
vorticity or pressure is negligible, the effects of the aether 
theory are indistinguishable from general relativity [6] . 
The aether theory can be tested only in high pressures 
and strong gravitational forces. The interior of neutron 
stars satisfies these conditions. In addition pulsars can be 
observed and studied empirically, enabling observational 
tests of theoretical models. 

Two sets of models define the structure of neutron 
stars. First models describing gravity that is the bind- 
ing force of the star, and second models describing the 
elementary constituents at the core of neutron stars and 
their repulsive forces that work against gravity and pre- 
vent the neutron star from collapsing and forming a black 
hole. Apart from quantum degeneracy pressure, strong 
nuclear interactions are the main sources of pressure in- 
side neutron stars. Various nuclear models give different 
pressure-density relations (equations of state) for the in- 
terior of neutron stars. Much of the uncertainty in the 
study of neutron stars is due to the lack of knowledge of 
the correct equation of state (EOS). Having this equa- 
tion and using a description of gravity we can find the 
mass-radius relation of neutron stars in static equilib- 
rium. In other words, for each neutron star of a given 
radius, we can find the mass for which the repulsive and 
attractive forces cancel and ensure the hydrostatic equi- 
librium of the star. In the context of general relativ- 
ity (or the gravitational aether theory), in contrast with 
Newtonian physics, this mass-radius relation has a max- 
imum mass M max beyond which no neutron star would 
exist and only black holes could have higher masses [8]. 
Therefore observations of high-mass neutron stars have 
the potential to constrain some equations of state and 
rule out others. Such an observation was made recently 
[9j. A millisecond pulsar was observed and its mass (1.97 
± 0.04 M ) was calculated using the Shapiro delay of the 
pulsar light due to its companion, a half solar-mass white 
dwarf. The high mass of this pulsar provides a lower limit 
on the maximum mass of neutron stars and rules out a 
number of proposed equations of state [TO] . Another ob- 
servation was reported in [TT]. The authors presented 
evidence that the black widow pulsar, PSR B1957+20, 
has a high mass. Their best fit pulsar mass was ~ 2.40 
± 0.12 M . A number of assumptions in the theoretical 
modelling of the pulsar contributed to the uncertainty in 



this number. Considering different constraints, the au- 
thors inferred a lower limit to the pulsar mass of M > 
1.66 M Q . Future observations of neutron stars will put 
additional constraints on the EOS. 



In this paper, we calculate the mass-radius relation 
predicted by the gravitational aether theory for two well- 
known equations of state. The first EOS (hereafter de- 
noted FPS) was calculated by Friedman and Pandhari- 
pande and improved by Lorenz, Ravenhall and Pethick 
[12] [13] . This equation of state is based on variational 
calculations over a wide density range using a realistic 
nuclear hamiltonian that contains two- and three-nucleon 
interactions, and fits the nucleon-nucleon scattering, as 
well as nuclear matter data. The Skyrme model is used 
in the FPS equation of state. In this model, the ef- 
fective interaction has the spatial character of a two- 
body delta function plus derivatives. The second EOS 
was calculated by Akmal and Pandharipande (hereafter 
denoted AP3)[II]. Some improvements of this calcula- 
tion compared to FPS are the use of Greens function 
Monte Carlo (GFMC) methods in the variational theory 
and including two-pion exchange three-nucleon interac- 
tion and isospin symmetry breaking terms in the hamil- 
tonian. Ozel and Psaltis have shown in [T5] that the 
complete mass-radius relation of neutron stars can be re- 
produced to high accuracy for all proposed equations of 
state, when the pressure of the neutron star is specified 
at three fiducial densities beyond the nuclear saturation 
density of p ns ~ 2.7 x 10 14 g cm ~ 3 . As they have cal- 
culated these values of pressure for the FPS and AP3 
equations of state, we will use their method to reproduce 
these two equations of state for densities higher than po 
which is a parameter to be adjusted for each EOS. For 
densities below p (the outer layers of the neutron star) 
we will use the SLy (Skyrme Lyon) equation of state cal- 
culated by Douchin and Haensel in [TS]. This equation 
is based on the effective nuclear interaction SLy of the 
Skyrme type, which is useful in describing the properties 
of very neutron rich matter. 



The structure of this paper is as follows: In section |TT} 
we summarize the phenomenological implications of the 



gravitational aether theory studied in [BJ. In section III 



we derive the equations of stellar structure for the grav- 
itational aether theory and relate the mass predicted by 
the theory to the observable mass of neutron stars. In 
section [TV] we solve the equations of stellar structure for 
a simplistic polytropic equation of state and explain the 
numerical method used. In section [V] we solve the equa- 
tions for the realistic FPS and AP3 equations of state, 
find the mass-radius relation of neutron stars predicted 
by the gravitational aether theory and compare it to the 
prediction of general relativity. An equivalent descrip- 
tion of the problem in terms of a modified EOS will be 
described in section IVT1 Section IVIII will include a discus- 
sion and the summary of our results. 
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II. PHENOMENOLOGICAL IMPLICATIONS OF 
THE AETHER THEORY 



III. THE AETHER EQUATIONS OF STELLAR 
STRUCTURE 



In [6], it was shown that the deviations of the aether 
theory from general relativity can only be significant in 
situations with rclativistic pressure, or (potentially) rela- 
tivistic vorticity. Furthermore, the authors demonstrate 
that the theory is consistent with all the current precision 
tests of gravity and cosmological observations. Here, we 
summarize their results: 

They show that for a perfect fluid with linear equa- 
tion of state (p oc p), the solutions to the gravitational 
aether theory are identical to those of general relativity 
only with a renormalized gravitational constant. As the 
gravitational coupling is not a constant in the aether, the 
authors find that in the case of homogeneous FLRW cos- 
mology, radiation energy gravitates more strongly than 
non-relativistic matter. The aether theory implies that 
gravity should be 33% stronger in the cosmological radi- 
ation era than the predictions of general relativity. 

As the increase of the gravitational constant at around 
the T = 0(1) MeV epoch induces an earlier freezeout of 
the neutron to proton ratio because of a speed-up effect of 
the increased cosmic expansion, the abundance of 4 He in- 
creases sensitively, and the abundance of deuterium (D) 
increases mildly while the abundance of 7 Be decreases. 
Comparing the theoretical prediction with the observa- 
tional light element abundances, the authors found that 
every light element abundance agrees with the gravita- 
tional aether theory within la. They found notably that 
7 Li fits the data better in the gravitational aether than 
in the standard big bang nucleosynthesis (which over- 
predicts 7 Li abundance by 4-5cr |17j). The main discrep- 
ancy found was with deuterium abundance observed in 
quasar absorption lines. 

Interestingly, cosmological observations, such as the 
CMB (Cosmic Microwave Background) [H] and the Ly-a 
forest [IH] prefer the aether prediction of stronger gravity 
in the radiation era, which is often interpreted as a larger 
effective number of neutrinos. 

The authors in [5] also examined the implications for 
precision tests of gravity using the PPN (parametrized 
post-Newtonian) formalism |20j . and showed that the 
only PPN parameter that deviates from its general rel- 
ativistic value is £4, the anomalous coupling to pressure 
(=1/3 for the aether and for general relativity), that 
has never been tested experimentally. Finally, they ar- 
gued that current tests of Earth's gravitomagnetic ef- 
fect mildly prefer a co-rotation of aether with matter, 
although they are consistent with an irrotational aether 
at 2cr level. 

In the current paper, we study the impact of the 
anomalous coupling to pressure in the aether on the 
structure of neutron stars and therefore suggest a novel 
way of putting constraints on the value of (4 using mass 
measurements of neutron stars. We will return to this 
point in section |V| 



As we assume a spherically symmetric static star, the 
metric will take the form: 
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with the line element being ds 2 = g^dx^dx". With this 
metric and since the problem is static the energy momen- 
tum tensor will be: 



( -e(r) 

p(r) 

p(r) 

\ p(r) 



where e(r) and p(r) axe the energy density and pressure 
at radius r of the star. The modified Einstein equation 
has the form: 



(5) 



where T M „ is given by: 



% v = (4/3) [T^ - -TSg^ + V{U^U V + g^)] , (6) 

Ta — 3p( r ) ~~ e ( r ) is the trace of the energy momentum 
tensor and V and are the pressure and four- velocity of 
the gravitational aether. Having Tjf = g^^Tp Vl where we 
have the Einstein summation over index /3, and imposing 
spherical static conditions on the aether we will get: 
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where: 



e r 



e(r) +p(r), 



p(r) = (1/3) [e(r)+p(r)]+(4/3)P(r). 



(7) 



As Eq. [5] is similar to the Einstein equation only with 
different energy density e(r) and pressure p(r) given by 
Eq. [7J the equations of stellar structure will be the same 
as in general relativity only with these updated quanti- 
ties (e and p). Using the Ricci tensor components and 
writing the different spherical components of the modi- 
fied Einstein equation (Eq. [5]) we will have [5T] : 
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where the prime superscript denotes the derivative with 
respect to radius. We have omitted the R^ equation as 
it is identical to Ree because of the spherical symmetry. 
Rewriting these equations using Eq. [7] we get: 
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Given our metric, the equation of hydrostatic equilib- 
rium for p and e is [ST] : 



~B 



5 + p 



(14) 



In addition to this, the same equation holds for our up- 
dated e and p. This equation is not independent from 
the modified Einstein equations and can be derived from 
them: 



&_ 2p' (e' +p')/4 + V 

B e+p~ e+p+V 



(15) 



Given the suitable boundary conditions, Equations [TT] 
12] 14 and 15 along with an equation of state giving e(p) 
(the energy density of the star as a function of its pres- 
sure) are enough to find our unknowns: A(r), B(r), p(r), 
e(r) and V(r). In practice, this needs to be done numer- 
ically. 

According to [7] the pressure in the vacuum does not 
vanish and is comparable to the pressure associated with 
dark energy. This pressure will be negligible for the cal- 
culations of neutron star structure. Therefore we can 
assume that pressure and energy density both vanish out- 
side the neutron star and the metric becomes the familiar 
Schwarzschild metric for which at r > R: 



B(r) = A^ir) = 1 - 



2GatM(R) 



(16) 



where R is the radius of the neutron star and M is the 
observed mass of the star given by the aether theory: 



M(R) = 
Using Eq. [7] this gives 
M(R) = 



e{r)r 2 drdQ. (17) 



[e(r) +p(r)]r 2 drdn. (18) 



The purpose of this work is to find the M — R relation 
for neutron stars assuming different equations of state. 
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FIG. 1. The p(r) — r relation of a neutron star of radius 
R = 1.175 with a polytropic equation of state p = Kp T where 
F = 9/5 for general relativity (solid) and the aether theory 
(dashed). The pressure of the aether V(r) which is negative 
is shown as well (dashed-dotted) . The units have been chosen 
so that Gjv = 1, c = 1 and K = 1. 



IV. NUMERICAL SOLUTIONS FOR A 
POLYTROPIC EQUATION OF STATE 

We start by solving the equations of stellar structure 
for a simple EOS. The polyropic equation of state is a 
power-law relation between pressure and matter density: 



p = Kp l 



(19) 



The energy density is given as a function of matter 
density by: 



e(p) = P c2 + P 



" pjp'W 

n P' 2 



(20) 



where the second term is negligible for non-relativistic 
matter. For neutron stars this term needs to be taken 
into account. If we use Eq. [l9]to perform the integration 
in the second term we will get: 



e(p) = pc 2 + 



r- r 



(21) 



If we choose our units so that c = 1 and Gjv = 1 (these 
will fix our time and mass units given a length unit) , the 
energy density as a function of pressure will be: 



y K J r-i 



(22) 



The differential equations [TT1 [12] [Ml and [15] along with 



Eq. 22 need to be solved numerically. The boundary 
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conditions are the values of A(R), B(R), p(R), e(R) and 
■p(R) (R is the radius of the neutron star). We set the 
values of pressure, energy density and aether pressure 
equal to zero at R: 



e(R) = 0, 
p(R) = 0, 
V(R) = 0. 



(23) 



The reason why we set the pressure of the aether equal 
to zero at the boundary of the star is that if the pressure 
of the aether in the vacuum is very small at infinity it 
will remain small up to the boundary of the star. This 
can be understood by writing equation |15| in the vacuum: 



which gives: 



V 



y/l - 2GM/r' 



(24) 



(25) 



If Voo — >• 0, then V — >• at the boundary of the star. It is 
straight-forward to generalize this argument to dynam- 
ical situations, i.e. aether pressure vanishes in vacuum 
everywhere if it vanishes at large distances. Therefore, 
the aether does not affect the binary mass measurements, 
enabling us to compare our results with the current ob- 
servations. 

The values of A(R) and B(R) are given by: 



B(R) = A -1 (R) = 1 



2M(R) 
R 



(26) 



As we do not have the value of M(R) and finding it is 
the purpose of this integration, we will use a shooting 
method in which for a given radius R we solve the dif- 
ferential equations with different values of M(R) starting 
from R/2 to smaller values. For each value of M(R) solv- 
ing the equations gives the energy density and pressure 
as a function of radius. Using Eq.[l8]we can find the inte- 
grated mass M. The value of M(R) for which M(R) = M 
is the correct mass of the neutron star of radius R. For 
this value the functions A(r), B(r), p(r), e(r) and Vir) 
are well-behaved and correspond to the solutions of equa- 
tions [TT] [12] [l4j [15] and [22] For instance, the pressure 
p(r) given by general relativity and the aether theory 
is shown in Figure [l] for a neutron star of fixed radius. 
The integration method used is a fourth-order Rungc- 
Kutta. Figure [2] shows the M(R) - R relation for neutron 
stars with a polytropic equation of state p = Kp T with 
r = 9/5 for general relativity and the aether theory in 
units for which Gjy = 1, c = 1 and K — 1. We have cho- 
sen r = 9/5 as it is consistent with the constraints found 
by Ozel and Psaltis in [15] and used in [22] to study the 
structure of neutron stars in /(R) gravity theories with 
perturbative constraints. We see that the aether theory 
gives a smaller mass for a neutron star of a given radius. 
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FIG. 2. The M(R) - R relation of a neutron star with a 
polytropic equation of state p = Kp T where F = 9/5 for 
general relativity (solid) and the aether theory (dashed). The 
units have been chosen so that Gn = 1, c = 1 and K = 1. We 
see that the maximum observed mass predicted by the aether 
theory is less than the maximum mass predicted by general 
relativity. We also see that in the Newtonian limit (large R) 
M(R) oc R- 7 . 



This was expected as gravity is stronger in the gravita- 
tional aether theory in the relativistic regime, reflected 
in the Q = 4Gat/3 relation. Therefore a neutron star of 
a given radius needs less mass to sustain its hydrostatic 
equilibrium compared to general relativity. 

To understand the behaviour of M(R) in large radii 
(the Newtonian limit) we can look at the equation of 
hydrostatic equilibrium in this limit: 



p(r)g(r) 



dp(r) 
dr 



(27) 



where g is the gravitational force and p(r) ~ e(r) in the 
Newtonian limit (as we have set c = 1). As p — Kp T the 
above equation takes the form: 
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We can write this as: 
M 



dr = (T- l)p L ~ 2 dp. 



(28) 



(29) 



In large radii we can treat M as constant and integrate 
both sides to get the following approximation: 



which gives: 
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FIG. 3 

densities p < po 



p — p (grey) and p— e/c (black) relations: Formatter 



lO^-Vcnr 1 we use the SLy EOS [16]. For 



p > po, the solid curves show the p — p (grey) and p — e/c 2 
(black) relations based on the minimal representation of the 
AP3 equation of state using the polytropic parameters of [15] . 
The dotted curves are the same relations for the FPS equation 
of state. The equations have been smoothed to avoid discon- 
tinuities in the derivative of pressure as the aether theory is 
sensitive to these derivatives. 



As r = 9/5 here we will have M oc R 7 , which is the 
behaviour seen in large radii in Figure [2j 



V. NUMERICAL SOLUTIONS FOR REALISTIC 
EQUATIONS OF STATE 

For densities below a fiducial density po of the order 
of the nuclear saturation density p ns ~ 2.7 x 10 14 g/cm 3 , 
the equation of state is well described by the SLy EOS 
[16j . For densities higher than po, it is shown in |15] 
that it is sufficient to specify the pressure of the neutron 
star at three fiducial densities p\ ~ 1.85p ns , p2 = 2pi 
and P3 = 2p2 in order to reconstruct the mass-radius 
relations based on the AP3 and FPS equations of state. 
Figure [3] shows the SLy equation of state in p < p as 
well as the parametrized AP3 and FPS equations of state 
in p > po using the polytropic parameters found by Ozel 
and Psaltis and listed in table [i] (from [15]). 

If we use the smoothed version of these two equations 
of state (to avoid discontinuities in the derivative of pres- 
sure as the aether theory is sensitive to these derivatives) 
to numerically solve the stellar structure equations in the 
same way we did for a polytropic EOS in section |TV] we 
find the M(R) - R relations shown in Figure [4] Gen- 
eral relativity predicts a maximum neutron star mass 



FIG. 4. (Color online) The M(R) - R relation of neutron 
stars given by general relativity (solid) and the aether theory 
(dashed) based on the parametrized AP3 (black) and FPS 
(grey) equations of state. The two observed pulsar masses of 
Demorest et al. [9] and van Kerkwijk et al. [11] are shown in 
orange and green respectively. 



M for the AP3 EOS. Therefore the difference in the two 
parametrized equations of state (as seen by looking at the 
difference of the solid and dotted curves in Figure |3| re- 
sults in a 24% difference in the maximum neutron star 
mass predicted by general relativity. The aether theory 
gives a smaller maximum mass as expected from Section 
IV The aether theory predicts M max ~ 1.58 M Q for the 
FPS EOS and M max ~ 2.00 M Q for the AP3 EOS. 

If the FPS equation of state is refuted on the basis of 
the 1.97 ± 0.04 M Q neutron star observed by Demorest et 
al. [9] and if we assume the validity of the AP3 equation 
of state, then both the aether theory and general relativ- 
ity agree with this observational measurement. For this 
EOS the maximum mass predicted by the Aether theory 
is 16% less than the maximum mass predicted by general 
relativity. However, for the AP3 EOS the aether predic- 
tion is inconsistent with the van Kerkwijk et al. [TT] mea- 
sured mass, although this mass measurement is uncertain 
due to assumptions made in the theoretical model used 
in calculating the neutron star mass. As other EOS can- 
didates exist (see [35]) more measurements of the mass 



TABLE I. Polyropic parameters for the FPS and AP3 equa- 
tions of state calculated by Ozel and Psaltis in [15] . 



Mn 



1.81 M for the FPS EOS and M n 



2.37 



EOS 


log po 


log Pi 


logP 2 


log P 3 


FPS 


14.30 


34.283 


35.142 


35.925 


AP3 


14.30 


34.392 


35.464 


36.452 
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FIG. 5. (Color online) The AP3 (solid black) and FPS (solid 
red) parametric equations of state p — e/c 2 and the modified 
effective equations of state p — e/c 2 given by the gravitational 
aether theory based on the same equations of state (dashed 
black for AP3 and dashed red for FPS). The blue line corre- 
sponds to p = e. 



If we use Eq. [7] to write e as a function of e and p and 
if we replace the derivatives with respect to radius with 
derivatives with respect to e we get: 



dp 
de 



P(£J 



rP-/(e) = 0, 



(33) 



where /(e) = dp(e)/de. Both /(e) and p(e) are given by 
the equation of state we are using. If we solve the dif- 
ferential equation [33] for the FPS and AP3 equations of 
state numerically we get p as a function of e. Using Eq. 
[7] we can find p as a function of e. These are the dashed 
curves shown in Figure [5] The figure compares the AP3 
and FPS parametric equations of state (p — e) and the 
modified equations of state (p — e) given by the gravita- 
tional aether theory based on the same equations of state. 
The effects of the aether theory become distinguishable 
beyond the nuclear saturation density (the region shown 
in Figure [5]). As in high densities, the aether theory gives 
a lower pressure than the one given by the equation of 
state, the stability of the neutron star will be obtained 
for a lower mass compared to the prediction of general 
relativity for a neutron star of the same radius. This is 
why the aether mass-radius relation falls below the gen- 
eral relativistic mass-radius relation for neutron stars. 



and radius of neutron stars are needed to put further con- 
straints on the existing equations of state, and allow us 
to make definite comparisons between the gravitational 
aether theory and general relativity. 

As was mentioned in section [TTJ, the value of £4 for 
the aether theory is 1/3 in contrast with general rela- 
tivity (C4 = 0). Therefore, assuming the AP3 equation 
of state, the maximum neutron star masses given by the 
aether theory and general relativity in Figure [3] can be 
translated into constraints on the value of £4. Using a 
linear interpolation, we find that £4 < 0.43 (0.19) at 95% 
confidence from the Demorest et al. [5] (van Kerkwijk et 
al. |llj ) mass measurements. 

VI. THE AETHER EQUATION OF STATE 

An equivalent description of the problem is suggested 
by Eq. [7] This equation gives the updated energy den- 
sity and pressure for which (with the new gravitational 
constant) the aether theory's Einstein equation ^ looks 
like the Einstein equation of general relativity. There- 
fore we can describe the aether theory's prediction of the 
structure of the neutron star as equivalent to the one of 
general relativity only with an updated equation of state 
given by p(e). To find this updated equation of state we 
equate the right-hand sides of equations fl4] and [TB"] t o get: 

_J^ = _^1 (32) 
e + p e + p 



VII. CONCLUSIONS AND FUTURE 
PROSPECTS 

The gravitational aether theory provides a possible so- 
lution to the cosmological constant problem. The struc- 
ture of neutron stars is related to their nuclear properties 
(the equation of state) as well as the theory describing 
gravity. Therefore, the aether theory can be tested in the 
light of mass-radius measurements of these stars. In this 
paper, we solved the aether theory's equations of stel- 
lar structure for two equations of state of nuclear matter 
(AP3 and FPS) and found the mass-radius relation of 
neutron stars based on these two EOS, and compared 
this with the mass-radius relation given by general rela- 
tivity. The FPS equation of state gave mass-radius rela- 
tions that both for general relativity and the aether the- 
ory, were incompatible with the 1.97 ± 0.04 M Q neutron 
star observed by Demorest et al. [5]. The mass-radius re- 
lations given by the aether theory and general relativity 
on the basis of the AP3 equation of state were both com- 
patible with this observation. We saw that for this equa- 
tion of state the aether predicts a maximum mass that is 
16% less than the maximum mass predicted by general 
relativity. We also found the modified equation of state 
of neutron stars given by the aether theory and based on 
that explained why the mass-radius relation given by the 
aether theory falls below the one given by general relativ- 
ity. It is important to note that there are other equations 
of state such as the one calculated by Miiller et al. [24] 
(and others mentioned in [53]) that do also agree with 
the Demorest et al. measured pulsar mass. 



In addition, including the effect of hyperons and quarks 
in the equation of state (e.g. [251430J ) can have a similar 
effect to the aether in lowering the maximum mass of neu- 
tron stars. For example, it is shown in [25] that a hybrid 
(nuclear+quark matter) star can have a mass-radius rela- 
tionship very similar to that predicted for a star made of 
purely nucleonic matter. The authors obtain hybrid stars 
as heavy as 2 M Q for reasonable values of their model pa- 
rameters. Due to these uncertainties in the equation of 
state we can not make definite comparisons between the 
aether theory and general relativity at the moment. 

To be able to test the aether theory more robustly, we 
need further constraints on the neutron star equations of 
state. In addition to constraints coming from the mass 
measurements of neutron stars (such as [S]), we also need 
further constraints from radius measurements that are 
considerably harder to get. The radius of a neutron star 
can be measured in various ways such as the thermal ob- 
servations of the surface of the neutron star, pulse profiles 
or redshift measurements |S]. We also know that merg- 
ers of compact objects such as pairs of neutron stars or 
neutron star-black hole pairs emit gravitational waves. 
These waves can be detected using current detectors if 
the emitter is close enough. Gravitational waves allow 
us to simultaneously measure masses and radii of these 
compact objects and could constrain the neutron star 
maximum mass and thus its equation of state [8] . They 
can also constrain the equation of state directly. In [3T] , 
the authors have studied the accuracy with which one 
can use gravitational wave observations of double neu- 
tron star inspirals to measure parameters of the neutron- 
star equation of state using numerical simulations. They 
concluded that gravitational wave observations could put 
a direct constraint on the EOS pressure at a rest mass 
density p = 5x 10 14 g cm~ 3 of 5p ~ 10 32 dyn cm~ 2 at 



an effective distance D c g = 100 Mpc (also see [32]). At 
this density, the difference of pressure between the aether 
EOS and the AP3 EOS is Ap = 7.6 x 10 32 dyn cm- 2 . 
This means that the aether theory's modified equation of 
state can be tested using gravitational waves, unless its 
predicted pressure value at the density mentioned above, 
is equal to the pressure predicted by another equation 
of state. This degeneracy will fade if further progress 
is made in constraining the equation of state of nuclear 
matter in densities above the nuclear saturation density. 
Another promising way to break this degeneracy is to 
study the dynamics of a collapsing neutron star which 
could distinguish the effects of modifying gravity from 
modifying the equation of state. 

Future observations ranging from the electromagnetic 
emissions of pulsars to the gravitational wave emissions of 
neutron stars in compact systems will allow us to learn 
not only about the nuclear constituents at the core of 
neutron stars but also about the nature of gravitation and 
fundamental questions such as the cosmological constant 
problem. 
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